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Abstract— We  consider  a  variation  of  the  Deffuant-Weisbuch 
model  introduced  by  Deffuant  et  al.  in  2000,  to  provide  new 
analytical  insights  on  the  opinion  dynamics  in  a  social  group. 
We  model  the  trust  that  may  exist  between  like-minded  agents 
through  a  trust  function,  which  is  a  discontinuous  (hard- 
interaction)  non-increasing  function  of  the  opinion  distance.  In 
this  model,  agents  exchange  their  opinions  with  their  neighbors 
and  move  their  opinions  closer  to  each  other  if  they  are  like- 
minded  (that  is,  the  distance  between  opinions  is  smaller  than 
a  threshold).  We  first  study  the  dynamics  of  opinion  formation 
under  random  interactions  with  a  fixed  rate  of  communication 
between  pairs  of  agents.  Our  goal  is  to  analyze  the  convergence 
properties  of  the  opinion  dynamics  and  explore  the  underlying 
characteristics  that  mark  the  phase  transition  from  opinion 
polarization  to  consensus.  Furthermore,  we  extend  the  hard- 
interaction  model  to  a  strategic  interaction  model  by  considering 
a  time-varying  rate  of  interaction.  In  this  model,  social  agents 
themselves  decide  the  time  and  energy  that  should  be  expended 
on  interacting  each  of  their  neighbors,  based  on  their  utility 
functions.  The  aim  is  to  understand  how  and  under  what 
conditions  clustering  patterns  emerge  in  opinion  space.  Extensive 
simulations  are  provided  to  validate  the  analytical  results  of  both 
the  hard-interaction  model  and  the  strategic  interaction  model. 
We  also  offer  evidence  that  suggests  the  validity  of  the  proposed 
model,  using  the  location  and  monthly  survey  data  collected  in 
the  Social  Evolution  experiment  over  a  period  of  nine  months. 
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1.  Introduction 

O  OBTAIN  a  fundamental  understanding  of  opinion 
formation  in  social  networks,  one  natural  question  to 
ask  is  that  how  the  initial  scattered  information  is  shared  and 
diffused  in  social  networks.  Indeed,  information  aggregation 
and  diffusion  generally  involve  interactions  between  agents 
in  the  network.  Given  the  role  of  social  interactions  in  the 
diffusion  of  information,  it  is  important  to  obtain  a  thorough 
understanding  of  how  the  structure  of  social  interactions 
affects  the  formation  of  opinion  and  in  shaping  individual 
behaviors. 

Indeed,  various  mathematical  models  of  social  interactions 
have  been  formulated  and  studied  historically  in  Economics 
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[l]-[3]  and  Social  Sciences  [4]-[7].  Studies  in  these  disci¬ 
plines  have  focused  on  providing  theoretical  tools  and  ana¬ 
lyzing  experimental  results  to  interpret  social  trends.  More 
recently,  an  increasing  number  of  disciplines  that  relate  to  in¬ 
formation  processing,  such  as  computer  science,  optimization 
and  control  as  well  as  signal  processing,  have  taken  a  close 
interest  in  the  study  of  social  interactions.  Another  motivation 
for  this  study  is  from  the  increasingly  dispersed  social  data 
in  clouds  of  interacting  computing  systems.  An  interesting 
question  is  whether  learning  and  parallel  decision  making 
in  such  environments  lead  to  herding  (a  phenomenon  that  is 
referred  to  as  information  cascade  [4],  [8])  or  fragmentation. 

A  popular  approach  in  Economics  and  Social  Sciences  is  to 
model  social  interactions  via  Bayesian  learning.  Specifically, 
Bayesian  models  describe  the  interactions  between  rational 
agents  by  postulating  that  agents’  actions  are  driven  by  the 
objective  of  maximizing  their  own  expected  utilities,  which 
depend  on  the  state  of  the  world  0.  Agents  observe  the  actions 
of  their  peers,  and  update  their  beliefs  using  Bayes  rule, 
optimally  fusing  public  information  obtained  by  observing 
their  peers’  behaviors,  and  private  information  about  0.  While 
it  is  possible  to  quantify  formation  and  propagation  of  opinions 
with  Bayesian  learning,  the  complexity  of  computing  the 
chain  of  actions  and  resulting  belief  updates  under  the  pure 
Bayesian  learning  model  complicates  the  analysis.  In  contrast, 
the  description  which  we  refer  to  as  non-Bayesian  learning, 
that  emerged  from  the  field  of  statistical  physics  [9],  radically 
reduces  the  complexity  of  the  interaction  model,  by  removing 
the  strategic  action  and  directly  postulating  a  rule  to  update 
agent  opinions.  In  this  case,  interactions  between  agents  are 
often  random  and  local  while  the  learning  rule  is  designed  to 
approximate  the  resulting  change  in  the  agents’  beliefs,  leaving 
out  the  agents’  decisions.  A  clear  exposition  of  the  difference 
between  the  two  approaches,  which  also  provides  an  excellent 
set  of  references  on  the  topics  is  in  [10]. 

The  interaction  models  discussed  in  this  paper  fall  in  this 
last  class  of  models  with  continuous  opinions  (see  [9]  for 
a  survey)  which  are  embodied  by  a  probability  distribution 
for  the  state  of  the  world  0  that  is  assumed  to  be  a  discrete 
random  variable.  Our  goal  is  to  present  a  simple  but  plausible 
model  for  the  evolution  of  opinion  in  a  population  of  agents, 
and  provide  answers  to  the  following  fundamental  questions 
concerning  the  formation  of  opinions  in  social  networks: 

•  What  kinds  of  interactions  will  lead  to  consensus  and 
which  ones  will  lead  to  polarization  of  opinions? 

•  Under  what  conditions  will  a  social  group  split  into 
clusters  of  opinions? 

In  the  following,  we  use  opinion  and  belief  interchangeably. 
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An  early  formulation  of  continuous  opinion  dynamics  was 
given  by  DeGroot  [11].  Indeed,  in  DeGroot’s  model,  individ¬ 
uals  start  with  an  initial  opinion  profile  on  a  binary  decision 
0  G  {0, 1},  represented  by  a  vector  of  the  probabilities  that 
each  agent  attributes  to  one  of  the  two  events,  say  Xi  =  P{0  = 
1),  which  is  agent  i’s  private  information.  The  update  process 
is  captured  by  a  fixed  stochastic  matrix  T,  where  the  (i, 
element  Tij  represents  the  relative  trust  that  agent  i  places  on 
agent  j’s  opinion.  Beliefs  of  individuals  are  updated  linearly 
by  taking  a  weighted  average  of  their  neighbors’  beliefs  (Tij 
are  the  weights).  The  analysis  of  such  a  model  is  relatively 
straightforward,  since  the  belief  updates  take  the  same  form  as 
the  well-known  average  consensus  algorithm  [12].  Therefore, 
considering  a  social  graph  Q  =  (V,  f )  with  V  representing  the 
set  of  agents,  and  such  that  £  includes  an  edge  {i,j)  if  and 
only  if  Tij  >  0,  consensus  of  opinions  is  attained  globally 
only  if  Q  is  connected,  otherwise,  each  connected  component 
will  attain,  in  general,  a  different  consensus  in  opinions. 

Another  related  model  is  the  Hegselmann- Krause  (HK) 
model  [13],  [14].  In  the  HK  model  as  well,  the  opinion  of 
each  agent  is  represented  by  a  real  number  0  <  Xi  <  1 
that  is  updated  synchronously,  similar  to  the  DeGroot  model. 
However,  the  HK  model  introduces  a  confidence  level  (or  a 
threshold)  tq  (set  to  1  in  [14])  to  model  the  lack  of  influence 
among  agents  whose  beliefs  are  too  far  apart.  Specifically, 
let  JVi  be  the  set  of  neighbors  (i.e.,  agents)  with  which 
agent  i  can  directly  communicate.  Let  Ii[k;To]  =  {j  G 
A/i  U i  :  \xi[k]  —  Xj  [k]  \  <  tq}  be  the  set  of  trusted  neighbors, 
i.e.,  whose  absolute  opinion  distances  are  less  than  tq  from 
agent  i’s  opinion  Xi[k]  after  the  interaction.  Then  after 
exchanging  beliefs  with  their  neighbors,  individual  opinions 
are  updated  synchronously  according  to  the  following  rule: 

'  jeii[k-,To] 

where  \S\  denotes  the  cardinality  of  the  set  S.  If  agent  i  has  no 
credible  neighbors  (i.e.,  Ii[k;  tq]  =  i),  then  it  does  not  change 
its  opinion.  Though  similar  to  the  DeGroot  model,  a  significant 
difference  is  that  the  HK  update  is  nonlinear  with  respect  to 
the  current  opinion  profile.  In  particular,  the  set  of  neighbors 
with  which  agent  i  updates  its  belief  may  change  with  k.  An 
extension  to  the  HK  model,  as  discussed  in  [13],  is  to  assume 
asymmetric  confidence  intervals  [— r^]  such  that  the  trusted 
set  li [k;  Tq]  is  replaced  by  li [k;  ri^Tr]  =  {j  G  A/i  Ui  :  —ri< 
Xj[k]  —  Xi[k]  <  Tr}.  Another  related  extension,  called  the 
heterogeneous  HK  (htHK)  model  [15],  is  to  introduce  diversity 
of  confidence  bounds  f  =  [ri,  •  •  •  ,  r^].  While  the  convergence 
properties  of  the  HK  model  can  be  observed  numerically, 
analytical  results  are  limited.  Recently,  [16]  extended  the 
analysis  of  the  HK  model  to  multidimensional  beliefs  and 
provided  an  upper  bound  on  the  convergence  time. 

The  Dejfuant-Weisbuch  (DW)  model  [17],  [18],  on  which 
our  proposed  analysis  is  based,  explores  the  effects  of  simple 
random  pair-wise  interactions  between  agents  whose  opinion 
distance  is  smaller  than  a  threshold.  In  this  model  as  well,  each 
individual’s  opinion  is  represented  by  a  real  number  Xj  G 
[0, 1].  Agents  i  and  j  are  randomly  selected  for  interaction, 
which  is  assumed  to  be  symmetric  (i.e.,  if  7  G  Mj,  then  j  G 
Mi).  Let  /i[/c;To]  =  {j  G  Mi  :  \xi[k]  -  Xj[k]\  <  tq}.  If 


j  G  Ii[k]TQ]  and  thus  i  G  /j[/c;ro],  then  after  the  interaction, 
opinions  are  updated  pair-wise  as  follows 

Xi  [k  +  l]  =  (1  -  fl)xi  [k]  +  flXj  [k]  (1) 

Xj[k  T  1]  =  (1  -  jTjXj [k]  +  flXi [k] , 

where  G  (0,  0.5]  is  called  the  mixing  parameter.  Specifically, 
Deffuant  et  al.  in  [17]  explored  this  system  over  a  square 
grid  in  which  individuals  are  only  connected  with  their  four 
immediate  neighbors.  Weisbuch  in  [18]  extended  this  simple 
lattice  topology  to  a  scale  free  network  topology.  Though 
different  from  the  HK  model,  the  DW  model  also  relies  on 
the  idea  of  bounded  confidence.  Both  models  focus  on  a  one 
dimensional  belief  so  they  exhibit  similar  behavior. 

The  model  proposed  in  this  paper  is  based  on  the  DW 
model,  but  rather  than  binary  decision  making,  as  assumed 
in  [13],  [14],  [17],  [18],  our  model  facilitates  the  analysis  of 
multi- alternative  decision  making  (decision  between  multiple 
alternatives),  similar  to  [16],  [19].  Specifically,  we  treat  each 
agent’s  opinion  as  a  vector  of  probabilities,  in  a  probability 
simplex  of  arbitrary  dimension;  each  element  of  the  opinion 
vector  represents  the  probability  that  a  certain  alternative  is 
true.  We  introduce  another  generalization  of  the  DW  model: 
a  state-dependent  trust  function  /i(d).  Although  similar  in 
spirit  to  the  parameter  /io  defined  in  [13],  [14],  [18],  [20], 

[21] ,  the  trust  function  jx{d)  in  our  model  varies  with  the 
squared  opinion  distance  d  between  the  interacting  agents. 
Clearly,  the  effect  of  //(d)  is  time  varying  since  agent  opinions 
evolve  over  time  and  its  value  depends  on  how  distance  is 
defined.  This  leads  to  the  so-called  hard-interaction  model 

[22]  that  includes  both  generalizations.  Moreover,  we  give  an 
explicit  mathematical  characterization  of  the  existence  of  a 
phase  transition  from  a  society  of  polarized  opinions  to  one 
with  a  convergent  opinion,  and  its  relation  to  the  society’s 
initial  opinion  profile.  Motivated  by  understanding  how  the 
rate  of  interaction  affects  the  outcome  of  the  interactions  in  the 
non-Bayesian  setting,  we  extend  the  hard  interaction  model  in 
[22]  by  considering  a  time-varying  rate  of  interaction,  which 
is  called  the  strategic  interaction  model.  In  this  model,  agents 
determine  how  much  effort  they  want  to  invest  in  interacting 
with  others.  Specifically,  at  each  time  instant,  agents  choose 
an  interaction  pattern  to  select  neighbors  (connected  agents), 
with  whom  interaction  produces  a  positive  net  profit.  Our 
focus  is  on  understanding  whether  simple  incentive  schemes 
for  interaction  can  lead  to  clustering  behaviors  in  opinion 
space.  Finally,  all  of  our  findings  are  validated  numerically. 
In  particular,  the  Social  Evolution  dataset  [23]  collected  from 
80  people  in  a  student  dormitory,  is  used  to  validate  the  model. 

H.  Social  Interaction  Model 

Let  V  =  -  ,n}  denote  a  set  of  social  agents,  who 

interact  with  each  other  over  a  fixed  undirected  communication 
graph  Q  =  (V,  f ),  where  £  is  the  set  of  edges.  We  focus  on  the 
case  where  Q  is  arbitrarily  connected  in  the  sense  that  there 
exists  at  least  one  path  connecting  any  two  agents  in  Q.  Denote 
Mi  the  set  of  agents  (also  called  neighbors)  connected  to  agent 
i  in  Q,  i.e.,  A/i  =  {j  G  V\i  |  (/,  j)  G  £}.  Agent  interactions  are 
modeled  as  pair-wise  random  encounters.  We  define  a  positive 
and  time-invariant  vector  p  whose  ith  element  pi  is  the 
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probability  of  agent  i  initiating  an  interaction  and  a  stochastic 
matrix  P  whose  (i,  element  Pij  is  the  probability  that 
agent  i  chooses  to  interact  with  agent  j.  We  assume: 

Assumption  1:  The  stochastic  matrix  P  is  time-invariant 
and  it  has  the  same  sparsity  as  the  undirected  network  graph 
Q,  i.e.,  Pij  >  0  VA  {i,j)  G  S. 

We  postpone  a  discussion  of  the  effect  of  time-varying  com¬ 
munication  rates  to  Section  IV. 

A.  Opinion  Updates 

The  opinion  evolution  runs  as  follows.  Each  agent  in 
the  network  starts  with  an  initial  opinion  profile  which  is 
modeled  as  a  g^-dimensional  vector  of  probabilities  Xj[0]  = 
[xji  [0],  •  •  •  ,  ^jg[0]]  in  opinion  space 

T'  =  <  X  =  [xi,  •  •  •  ,  Xq]^  I  =  1  and  xi  G  [0, 1]  >  . 

I  £=1  J 

The  opinion  Xj  [0]  can  be  thought  of  as  a  set  of  probabilities 
that  the  agent  attributes  to  outcomes  of  a  discrete  random 
variable,  or  a  mixed  strategy.  At  each  step,  the  network  selects 
an  agent \  say  agent  i,  with  probability  pi  to  initiate  an 
interaction;  and  with  probability  Pij,  agent  i  will  then  choose 
to  interact  with  agent  j  ^  Mi.  If  neither  of  the  agents  finds  the 
opinion  of  the  other  agent  trust-worthy  (according  to  a  trust 
function),  then  nothing  happens.  Otherwise,  the  agent  who 
is  open-minded  to  the  other  agent’s  opinion,  will  update  its 
opinion  in  a  way  that  decreases  the  opinion  distance  between 
the  two  agents.  Motivated  by  the  application  in  [24],  we 
introduce  a  distance  function  s(xi,Xj)  :  X  x  X  — ^  M+  to 
measure  the  difference  between  opinions,  i.e., 

s(Xj,Xj)  =  ||Xj  -XjlU  =  [(Xi  -Xj)^^(Xi  -XjU/2 

where  A  G  is  a  positive  definite  matrix.  The  set  X  is 

bounded  with  respect  to  the  norm  ||  x,|U  :=  Ax.,0),  i.e., 
supj|xi||^  <  oo  for  Mxi  G  A'  and  Vi  G  V.  Hence,  the 
triangular  inequality  implies  s{xi,Xj)  <  2supJ|x^m  :=  Sgup, 
where  sup^  denotes  the  supremum  of  the  value  x^  for  Vi  G  V. 

For  ease  of  notation,  let  dij[k]  =  s‘^{xi[k],Xj[k])  denote 
the  squared  opinion  distance  after  k  network-wide  interactions 
have  occurred  and  dgup  =  denotes  the  supremum  of 
squared  distance.  Given  that  agents  i  and  j  interact  at  the 
{k  step,  the  opinions  are  updated  as  follows: 

Xi  [kpl]=  Xi  [k]  P  /i  {dij  [k] )  {xj  [k]  -  Xi  [k] )  (2) 

Xj[kpl]  =xij[k]  P  p{dij[k]){xi[k]  -  xij[k])  (3) 

where  p{d)  defines  the  trust  function  of  each  agent  in  terms 
of  d.  Note  that  //(d)  >  0  when  the  agent  is  open-minded  to 
the  opinion  of  the  interacting  agent  and  p{d)  =0  otherwise. 

We  say  that  the  network  attains  consensus  (herding)  when 
all  agents  have  the  same  opinion  vector  with  respect  to  the 
distance  measure,  i.e.,  \fi,j  G  V,  s(xi,Xj)  =  =  0, 

but  need  not  believe  in  only  one  outcome,  i.e.,  consensus 

^The  phrase  ‘randomly  selected  by  the  network’  is  not  meant  to  imply 
that  the  ’’network”  chooses  the  initiating  agent;  rather  this  is  a  model  for 
agents  randomly  getting  activated.  In  earlier  work,  we  had  used  the  model  of 
a  common  rate  Poisson  clock  at  each  agent,  which  dictates  when  the  agent 
initiates  an  interaction 


does  not  require  Xie[k]  =  1  for  \/i  and  some  0  G  [1,. . .  ,q]. 
In  contrast,  polarization  describes  a  process  in  which  the 
society  is  asymptotically  divided  into  non-interacting  sub¬ 
groups  (herds),  each  of  which  is  internally  in  consensus  with 
respect  to  the  distance  measure.  More  precisely,  s(xi,Xj)  =  0 
if  i,j  belong  to  the  same  sub-group.  Finally,  we  say  that  a 
system  devolves  into  clusters  of  opinions  if  there  exist  distinct 
groups  of  agents  with  similar  opinions  within  the  group,  while 
the  opinion  distances  between  different  clusters  are  relatively 
large.  Hence,  the  opinion  distances  between  agents  in  the  same 
cluster  need  not  be  zero.  Note  that  herding  or  clustering  is  in 
opinion  space;  it  does  not  make  any  assumption  about  Q. 

The  dynamics  of  individual  opinions  clearly  depend  upon 
li(d).  We  present  here  one  specific  construction  of  //(d). 

B.  Hard-Interaction  Model 

The  hard-interaction  model  considers  a  trust  function  that 
satisfies  the  following  assumption. 

Assumption  2:  The  trust  function  p{d)  G  [0,0.5]  is  a  non¬ 
increasing  function  of  the  squared  opinion  distance  d. 

The  function  p(d)  models  the  trust  that  may  exist  between 
like-minded  agents.  Given  that  the  amount  of  trust  one  places 
in  another  reflects  the  amount  of  risk  one  is  willing  to  take 
in  social  interactions,  it  is  reasonable  that  as  the  difference 
between  the  opinions  of  two  interacting  agents  increases,  the 
trust  p{d)  between  the  agents  decreases.  Simply  put,  agents 
interact  more  strongly  with  agents  that  have  closer  opinions. 
Note  also  that  p(d)  is  similar  in  spirit  to  the  fixed  mixing 
parameter  fl  defined  in  [13],  [14],  [18],  [20],  [21],  but  here 
it  is  dependent  on  the  opinions  of  the  two  agents.  Hence, 
p{d)  is  stochastic:  its  value  depends  on  which  two  agents  are 
communicating  and  on  their  opinion  distance.  The  constraint 
p(d)  <  0.5  implies  that  an  agent  trusts  its  own  prior  opinion 
at  least  as  much  as  that  of  the  opinion  of  its  neighbor.  The 
technical  reason  will  become  apparent  soon. 

The  hardness  of  the  hard-interaction  model  stems  from  the 
existence  of  a  threshold  r,  such  that  if  the  squared  opinion 
distance  between  interacting  agents  is  larger  than  r,  then 
neither  of  the  two  agents  updates  its  opinion. 

Assumption  3:  There  exists  a  threshold  r  :  d  >  r  ^ 
p(d)  =  0  and  p(0) / p(r~)  <  P  <  oo. 

This  implies  that  agents  will  interact  with  their  neighboring 
agents  if  and  only  if  d  is  less  than  a  given  r  G  (0,(isup). 
Therefore,  the  system  in  (2)  and  (3)  might  not  change  at  all 
after  an  opinion  exchange,  implying  that  the  agents  have  no 
influence  over  each  other  when  they  are  not  like-minded^. 
For  reasons  that  will  soon  become  apparent  (Femma  6),  even 
though  our  results  are  not  restricted  to  the  step-function  in  (1), 
the  condition  p(0)/p{r~)  <  oo  means  that  there  must  be  a 
discontinuity  of  p{d)  at  d  =  r  for  the  results  to  hold. 

Assumption  4:  The  trust  function  p(d)  is  concave  and  C^- 
differentiable  for  Vd  G  (0,  r). 

The  regularity  condition  imposed  in  Assumption  4  is  needed 
for  analytical  reasons,  made  clear  in  the  next  section. 

^The  system  will  not  change  if  agents  are  already  in 

agreement. 
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Fig.  1.  Opinions  of  agents  i^j  and  i  before  and  after  the  interaction  of 
agents  z,  j,  where  x'  and  s'  =  ^/d/  denote  the  opinion  and  opinion  distance 
after  a  generic  update. 

With  all  of  the  above  in  mind,  it  can  be  checked  that  the 
DW  update  in  (1)  is  a  special  case  of  our  proposed  model  in 
(2)  and  (3).  Let  f{x)  be  the  step  function.  It  follows  from  (1) 
that  the  squared  opinion  distance  after  the  update  is  given  by 

\xi[k 1]  —  Xj[k 1]\‘^  (4) 

=  [1  -  2/7/  -  \xi[k]  -  Xj[k]\^)]‘^  ■  \xi[k]  -  Xj[k]\^. 

Let  Xi  =  [xi,  1  —  Xi]^  for  Vi  G  V.  Then  with  q  =  2,  A  =  /2/2, 
fi{d)  =  /i  for  0  <  d  <  Tq  ,  our  model  in  (2)  and  (3)  reduces 
to  the  DW  model  in  (1)  in  terms  of  the  change  in  squared 
opinion  distance.  Note  also  that  for  the  DW  case,  the  concavity 
condition  (Assumption  4)  is  met. 

III.  Opinion  Dynamics 

In  this  section,  we  study  the  asymptotic  convergence  prop¬ 
erty  of  the  opinion  dynamics  under  the  modeling  assumptions 
described  in  the  previous  section.  Our  aim  is  to  provide 
analytical  insights  on  what  may  happen  as  the  initial  scattered 
information  is  diffused  in  the  network,  with  an  emphasis  on 
finding  global  conditions  under  which  a  society  will  converge 
to  one  opinion  or  it  will  polarize  into  isolated  groups. 

Suppose  that  the  (/c  -f  1)^^  interaction  is  between  agents  i 
and  j.  It  follows  from  (2)  and  (3)  that 

dij[k  +  1]  =  [1  -  2^i{dij[k])f  dij[k]  (5) 

=  dij[k]  -  4:ij{dij[k])  (1  -  iJ,{dij[k]))  dij[k]. 

Thus  the  condition  |1  —  2fi{dij)\  <  1  (see  Assumption  2) 
is  needed  to  ensure  that  distances  do  not  increase  after  an 
interaction.  We  call  this  decrease  in  dij  the  private  marginal 
benefit  that  is  caused  by  the  interaction  between  agents  i  and 
j.  Note  that  the  condition  1/2  <  p{dij)  <  1  also  leads  to 
a  decrease  in  dij,  but  the  “ordering”  of  the  opinions  and 
xj  reverses  after  the  interaction,  implying  that  the  agent  puts 
more  trust  in  the  prior  opinion  of  its  neighbor  than  in  its  own. 

Because  of  the  triangular  inequality  (see  Fig.  1),  if  opinions 
of  the  interacting  agents  {i,j)  are  moving  closer,  then  the 
perimeter  of  the  triangle  with  vertices  i,j  and  any  other  point 
i  shrinks  after  the  update^.  In  other  words,  the  sum  of  opinion 
distances  can  only  shrink  after  an  interaction;  any  interaction 
between  the  pair  (i,  j)  G  E  also  changes  the  opinion  distances 
between  these  two  nodes  and  their  neighbors.  This  creates  a 
network  effect  called  network  externality. 

Now  consider  the  dynamic  effects  that  arise  from  network 
externality.  Since  opinions  of  the  non-interacting  agents  do 

^Of  course,  this  does  not  imply  that  the  opinion  distances  of  ^  to  z  and 
j  both  shrink.  It  is  quite  possible  that  one  (but  not  both)  of  the  interacting 
agents  moves  away  from  i. 


not  change  after  the  interaction,  i.e.,  Xi[k  +  1]  =  Xi^[k]  for 
Ml  ^  i,j,  di  simple  manipulation  of  (2)  and  (3)  yields 

du[k  +  1]  =  du[k]  +  p  {dij[k])dij[k] 

+2^i{dij[k]){-x.i[k]  -  -x.i[k\f  A{yLj[k]  -  Xi[k])  (6) 

dje[k  +  1]  =  dje[k]  +  ij,‘^{dij[k])dij[k] 

+2ii{dij[k]){xj[k]  -  xe[k])'^ A{xi[k]  -  Xj[k]).  (7) 

Let  D^^k  +  1]  :=  Y17=iY^m=r+i<^rm[k  +  1]  be  the  sum 
of  squared  distances  over  all  possible  pairs  of  agents  in 
the  network,  given  that  agents  i  and  j  interacted  at  /c  +  1. 
Let  D[k]  =  Ylr=iYlm=r+i^rm[k]  be  the  sum  of  squared 
distances  at  time  k  which  is  prior  to  the  interaction.  Thus 
D[k]  is  a  measure  of  the  disparity  of  opinion  in  the  society. 
The  change  in  D  after  the  interaction  equals  (see  Appendix) 

D^^[k  +  1]  -  D[k]  =  -2np{dij[k])  [1  -  p{dij[k])]  dij[k]. 

(8) 

Notice  that  the  change  in  the  overall  sum  of  squared  distances 
depends  entirely  on  dij  [k]  between  the  two  interacting  agents 
and  it  is  n/2  times  as  large  as  the  change  in  dij[k]  of  the 
interacting  pair  (see  Eqn.  (5)). 

A.  ODE  Approximation  of  the  Distance  Dynamic 

Agents  in  the  network  interact  at  random:  at  each  time, 
agent  i  is  selected  with  nonzero  probability  pi  to  initiate  an 
interaction  and  with  probability  Pij,  agent  i  will  then  choose 
to  interact  with  agent  j.  Let  P  G  be  the  matrix  of 

probabilities  of  the  pair  {i,j)  performing  any  exchange,  i.e., 
Pij  =  Pi  Pij  PpjPji.  The  first  term  corresponds  to  the  case  of 
agent  i  initiating  an  interaction  with  agent  j  and  the  second 
to  that  the  probability  of  agent  j  initiating  the  interaction  with 
agent  i.  By  Assumption  1,  P  =  [Pij]  has  the  same  sparsity 
as  Q,  i.e.,  Pij  >  0  ga  {i,j)  G  £.  Hence,  given  that  Q  is 
connected,  information  can  be  propagated  directly  or  indirectly 
from  any  agent  Q  to  any  other  agent  with  nonzero  probability. 

Let  fk  {d)  be  the  distribution  of  d  between  any  pair  of  agents 
at  time  k.  For  n  sufficiently  large,  the  conditional  expectation 
of  d[k  +  1]  with  respect  to  fk-\-i{d)  can  be  approximated 
by  the  sample  mean  and  we  have  the  following  relation: 
EMfc  +  l]}  =  E{d[k  +  1]|  I  {i,j)  interacts}  ^ 

^(^ij)^s  -Pzj  rf( J-i)/2  •  With  this  sample  mean  approximation 
and  the  relation  in  (8),  we  have 

E{d[kpl]}  -E{d[k]} 

_  Di^[k  +  l]-D[k] 

^  n(n  — l)/2 

=  —  ^  ^ij E{dij[k])  [1  —  k'idijlk])]  dij[k]. 

Our  general  approach  to  studying  the  asymptotic  behavior 
of  opinion  dynamics  here  is  to  map  the  difference  equation 
onto  an  ordinary  differential  equation  using  an  Euler-type 
approximation,  i.e.,  the  derivative  of  a  continuous  function 
d{t)  equals  d{t)  ^  {d[k  +  1]  —  d[k])/h  where  h  >  0  denotes 
the  discretization  step  size  and  d[k]  is  the  value  of  d{t)  at 
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t  =  kh.  Using  Euler’s  approximation  and  setting  h  =  the 
following  ordinary  differential  equation  (ODE)  can  be  derived 


^^l 

o 


d{t)  =  -  Aj  Kdjj  (0)  [1  -  {t))]  dij  (t)  (9) 

=  —  'y  ^  Pijp{dij(t))dij(t). 

where  d{t)  =  K{d{t)}  denotes  the  mean  of  the  average 
distribution  ft{d).  Eor  convenience,  we  do  not  explicitly  show 
the  time  variable  t  whenever  it  does  not  cause  confusion  in 
the  rest  of  the  paper. 

Taking  the  expectation  on  both  sides  of  Eqn.  (9)  with  respect 
to  the  distribution  /t(d),  we  get 


<i=  -  y]  E 


Pij  p{dij  ^dij 


(10) 


E  \  pi^dij^dij  /  —  ^  {p{dij^dij^  . 


(ij)e£ 

The  second  equality  holds  true  because  it  follows  from  As¬ 
sumption  1  that  given  a  network  Q,  the  pairwise  communica¬ 
tion  rates  Pij  are  fixed  prior  to  any  interaction.  We  are  now 
ready  to  develop  the  main  results  regarding  opinion  dynamics 
in  our  model. 


B.  System  Analysis  and  Results 

The  next  property  gives  a  lower  bound  on  the  dynamic  of 
the  expected  squared  distance  d. 

Lemma  5:  Under  Assumptions  2-4,  the  dynamic  in  (10) 
is  lower  bounded  by  d  >  —f3p{d)d  . 

Proof:  See  Appendix.  ■ 

It  is  to  be  noted  that  Assumptions  2-  4  allow  us  to  use  the 
expected  squared  distance  d{f)  as  a  surrogate  for  studying  the 
evolution  of  total  disparity  D{t)  between  agents’  opinions. 
Erom  Lemma  5,  it  is  expected  that  (10)  will  not  converge  if  the 
lower  bound  does  not  converge.  Eor  convenience,  we  express 
the  dynamic  of  the  lower  bound  system  by  b  =  —I3p{b)b. 
Indeed,  the  ODE  of  b  locally  resembles  the  form  of  the  logistic 
equation,  as  initially  investigated  in  [22] .  However,  the  model 
in  this  paper  includes  network  externalities  that  were  not 
accounted  for  in  [22]  and  hence  analysis  is  quite  different. 

Lemma  6:  Under  Assumptions  2  and  4,  the  system  b  = 
—Pp{b)b  converges  if  and  only  if  r  >  6(0). 

Proof:  See  Appendix.  ■ 

Using  Lemma  6,  a  necessary  condition  for  (10)  to  converge 
under  the  hard-interaction  model  can  be  established  and  it 
generalizes  the  result  in  [22]  to  the  model  with  externalities 
as  defined  in  this  paper. 

Lemma  7:  Under  Assumptions  1-4,  when  the  number  of 
agents  n  is  sufficiently  large,  a  necessary  condition  for  the 
interaction  model  in  (2)  and  (3)  in  a  connected  network  to 
converge  almost  surely  is  r  >  d[0]. 

Lemma  7  indicates  that  if  the  system  converges  almost 
surely,  then  the  threshold  must  be  above  d[0],  the  average 
initial  squared  distance.  If  r  <  d[0],  then  agents  will  remain 


Fig.  2.  Group  polarization 


polarized  with  a  positive  probability.  One  way  to  interpret 
the  threshold  value  r  is  that  it  represents  how  open-minded 
a  society  is.  If  it  is  open-minded  enough  (relative  to  this 
initial  dissonance  in  opinions),  then  the  system  will  converge. 
This  also  implies  the  existence  of  a  phase  transition  Tpt 
from  polarization  to  consensus  for  Tpt  sufficiently  above  d[0]. 
However,  we  will  show  next  that  the  necessary  condition 
developed  in  Lemma  7  is  not  a  sufficient  condition. 

Let  there  be  two  groups  of  opinions  Hi  and  H2  in  a  society 
(See  Eig.  2).  Within  each  group,  the  opinions  are  in  consensus, 
but  for  V(i,  j)  :  i  G  i7i,  j  G  H2,  the  distance  dij  =  r  +  e  is 
above  the  threshold.  Let  £12  =  {(^,7)  |  i  G  i7i  and  j  G  H2}. 
Since  the  two  groups  do  not  communicate  (as  the  trust  function 
p{dij)  =  0  :  i  G  i7i,7  G  H2),  the  agents  cannot 

come  to  a  consensus.  However,  if  the  number  of  intergroup 
pairs  of  agents  normalized  by  the  total  number  of  pairs  is 
n(i^'i)/2  <  tG’  then  d[0]  =  (r  +  An(i-i)/2  <  ^  and  thus 
p  (d[0])  >  0.  This  contradicts  the  assumption  that  the  two 
group  do  not  communicate.  Hence,  r  >  d[0]  is  not  a  sufficient 
condition  for  consensus. 

A  final  remark  pertains  to  the  effect  of  Assumption  1  on 
the  network  dynamics.  Recall  that  in  the  derivation  of  (10), 
we  have  shown  that  the  rate  of  interaction  Pij  does  not  affect 
the  dynamic  of  the  mean  d{t)  when  Pij  is  fixed  (as  required 
in  Assumption  1).  If  this  is  the  case,  then  for  n  sufficiently 
large,  there  is  no  graph  selection  that  can  mitigate  the  effect  of 
a  large  initial  average  distance  d[0]  (relative  to  the  threshold) 
since  the  necessary  requirement  for  consensus  in  Lemma  7 
does  not  change.  Then  the  natural  question  to  ask  is  that  what 
asymptotic  opinion  profile  can  be  observed  if  Pij  is  time- 
varying.  We  will  examine  this  in  the  following  section. 

IV.  Strategic  Opinion  Eormation 

While  the  interaction  model  discussed  in  Section  III  studies 
the  opinion  dynamics  under  random  interactions  with  a  fixed 
rate  of  communication,  there  are  many  social  settings  in 
which  agents  dynamically  determine  how  much  effort  and 
time  they  want  to  consume  in  interacting  with  their  neighbors. 
A  strategic  model  of  interaction  for  an  agent  is  the  process 
of  defining  a  goal  through  a  utility  function  and  allocating 
resources  (i.e.,  time  and  energy)  for  interactions  with  its 
neighbors  to  maximize  the  utility  function.  In  this  model,  we 
contemplate  the  possibility  of  agents  not  wanting  to  interact  at 
all  with  their  peers  and,  therefore,  we  allow  the  sparsity  pattern 
of  P  to  be  different  from  that  of  Q  and  the  communication 
rates  Pij  for  V(7,7)  G  f  to  vary  with  time. 

There  are  two  steps  in  modeling  strategic  interactions  for  a 
network  of  agents.  Eirst,  one  needs  to  explicitly  model  the 
cost  and  benefit  that  each  agent  receives  as  a  result  of  a 
particular  action,  which  define  the  incentives  agents  have  to 
interact  more  or  less  often  with  neighboring  agents.  In  this 
case,  it  is  assumed  that  the  initiator  of  an  interaction  knows 
the  opinions  of  all  its  neighbors,  before  it  starts  the  interaction. 
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Second,  the  strategic  model  should  be  tractable  so  that  it  can 
provide  insights  or  predictions  on  how  individual  incentives 
affect  the  aggregation  of  information  and  the  formation  of  the 
asymptotic  opinion  profile. 

A.  Strategic  Interaction  Model 

As  highlighted  earlier,  the  utility  function  of  an  agent  should 
reflects  its  goal  for  interactions.  The  mathematical  description 
of  goal-attainment  thus  involves  the  cost  and  benefit  for  par¬ 
ticipating  in  a  particular  interaction.  We  define  the  following 
utility  function:  for  some  k>  1, 

G ([AilieM )  =  E  Pii  ^ )]  - 

jeMi 

where  r{dij)  is  the  benefit  or  reward  agent  i  receives  from 
interacting  with  agent  j  and  ^{dij)  represents  the  cost  of 
interacting  with  agent  j.  Both  the  cost  and  reward  are  functions 
of  the  squared  opinion  distance.  Note  also  that  without  the 
second  term  ,  maximizing  the  utility  function  is  equivalent 
to  assigning  Pij  =  1  to  agent  j  that  gives  the  highest  positive 
profit,  measured  by  the  reward  received  minus  the  cost.  It  is 
however,  important  to  emphasize  that  the  focus  of  the  strategic 
interaction  is  not  to  model  the  agents  as  selfish  individuals  who 
calculate  their  overall  profits  from  each  potential  interaction, 
but  rather  to  study  their  tendency  to  initiate  interactions  that 
are  beneficial  to  them.  This  tendency  is  measured  by  Pij  and 
the  term  P-^  in  (11)  should  be  interpreted  as  the  force  behind 
such  incentives.  The  mathematical  formulation  of  the  utility 
function  thus  allows  us  to  translate  individual’s  tendency  to 
interact  into  a  tractable  process  that  can  potentially  provide 
answers  to  why  herding  behavior  occurs. 

Recall  that  the  underlying  assumption  is  that  agents  do  seek 
to  arrive  at  consensus,  but  subject  to  the  interaction  rules. 
Also  recall  from  Eqn.  (8)  that  the  private  marginal  benefit  (5) 
enjoyed  by  the  interacting  pair  (i,  j)  is  amplified  n/2  times  (n 
is  the  number  of  nodes)  in  terms  of  the  network-wide  average 
squared  distance.  Motivated  by  this,  we  set  the  reward  function 
to  be  proportional  to  this  private  marginal  benefit,  i.e., 

r{dij)  =  ap{dij)dij, 

where  p{d)  is  defined  in  (9)  and  a  is  called  the  reward 
coefficient. 

Let  us  now  consider  the  cost  of  interaction  f{dij),  which 
should  specify  the  energy  for  communicating  with  agent  j 
whose  opinion  is  away  from  agent  i.  Here  we  present 

two  interesting  constructions  of  the  energy  function 

(i)  f{dij)  =  ^idij  if  dij  <  T  and  f{dij)  =  +oo  if  dij  >  r 

(ii)  f{dij)  =  ^2P^{dij)dij  if  dij  <  r  and  f{dij)  =  +oo  if 

dij  >  T 

where  71  and  72  are  called  the  cost  coefficients  associated  with 
the  two  energy  functions.  The  first  construction  of  the  energy 
function  assumes  that,  if  the  (squared)  opinion  distance  dij  is 
large,  then  agent  i  consumes  more  energy  to  communicate  to 
j  than  when  dij  is  small.  The  notion  here  is  that  it  will  take 
more  effort  to  dialogue  with  someone  farther  away  in  opinion 
and  convince  them  to  move  their  opinion  profile  closer.  Rather 
than  choosing  a  linear  cost  function,  we  could  have  chosen 
a  function  of  the  form  f{dij)  =  ^if{dij)dij  where  f{d)  is 


any  positive  non-decreasing  function  of  d.  As  will  become 
evident  from  the  analysis  in  the  next  section,  the  choice 
f{d)  =  l  does  not  lead  to  any  loss  of  generality.  In  contrast, 
the  second  construction  of  the  energy  function  captures  the 
overall  cost  of  opinion  changes  for  both  agents.  Note  that 
given  dij  <  r,  the  amount  of  energy  spent  in  interacting  with 
agent  j  is  proportional  to  agent  j’s  change  in  opinion,  i.e., 
d{xj  [k] ,  Xj  [k-\-l])  =  jjf  {dij  [k])dij  [k] ,  which  also  equals  agent 
Es  change  in  opinion  d{xi [k] ,  Xi[k 1]) .  In  both  cases,  agents 
assign  infinite  energy  for  interactions  with  neighbors  who  are 
not  like-minded.  Each  of  the  two  cases  will  lead  to  different 
opinion  formation  processes,  as  we  shall  see  shortly. 

Agents  interact  with  each  other  according  to  the  following 
strategic  model.  An  agent  is  randomly  activated  at  each  time 
instant  with  equal  probability,  i.e.,  pi  =  1/n  for  Vi  G  V.  If 
agent  i  is  selected,  it  needs  to  make  a  strategic  decision  on 
an  interaction  pattern  for  selecting  neighbors.  The  decision  is 
in  the  form  of  probabilities  Pij  obtained  by  maximizing  the 
utility  function  defined  in  (11).  Then  based  on  the  decision, 
agent  i  interacts  with  one  of  the  neighbors  in  Mi,  and  this 
event  is  followed  by  the  opinion  update  rule  given  in  (2)  and 
(3).  At  the  next  time  instant,  the  same  procedure  is  applied. 


B.  System  Analysis 

In  this  section,  we  explore  how  the  individual  incentives  for 
interactions  affect  the  opinion  formation  process.  Suppose  that 
agent  i  is  chosen  at  time  k.  Recall  that  the  first  step  in  the 
strategic  interaction  model  is  to  determine  Pij  for  Vj  G  Mi 
maximizing  the  utility  function,  i.e., 

max  E  Pij  -  ^{dij[k])]  - 

under  the  constraint  that  ^j^j^.^i  Pij  =  1  and  Pij  >  0. 
Let  [a]+  =  a  if  a  >  0,  and  0  otherwise.  Solving  the  above 
optimization  with  respect  to  Pij  yields,  for  Vj  G  Mi 

Pii  [k]  =  J-  {Wii  [k])  -  adij  [fc])]+)  “  (12) 

where  the  scaling  factor  Si  is  defined  as  follows:  let  c  = 

{[r{di^[k])  -  adiM]^)~\ 


Si  := 


if  c  7^  0 
otherwise. 


(13) 


In  the  case  that  Pij[k]  =  0  for  Vj  G  Mi,  we  have  Pii[k]  =  1 
and  Pii[k]  =  0  otherwise.  That  is,  if  interacting  with  any  of 
its  neighbors  will  bring  zero  or  negative  net  benefit  to  agent 
i,  then  agent  i  will  choose  to  be  a  recluse  for  the  moment. 

1)  Case  Study  (i):  Consider  the  case  when  the  energy 
function  is  defined  as  f{d)  =  yid  for  d  <  r  and  +(X) 
otherwise.  Note  that  under  Assumption  2,  the  term  p{d)  = 
/i(d)(l  —  ii{d))  is  a  non-increasing  function  of  d,  with  maxi¬ 
mum  value  p(0).  When  dij  <  r,  replacing  f{dij)  with  yidij  in 
(12),  we  observe  that,  if  the  ratio  yi/a  of  the  cost  coefficient 
to  the  reward  coefficient  is  greater  than  or  equal  to  the  product 
p(0),  then 


(14) 
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for  ^dij .  This  situation  implies  that  the  agents  have  insufficient 
incentives  to  interact  with  each  other.  In  this  case,  it  follows 
from  equation  (12)  that  Pij  =  0  for  Vj  G  A/^,  given  that 
dij  <  r.  Thus  agents  always  choose  to  not  interact  with  the 
other  agents  and  remain  inactive,  i.e.,  Pii[k]  =  1  for  V/c.  As  a 
result  from  (14)  and  the  fact  that  Pij  =  0  whenever  dij  >  r, 
the  system  will  not  evolve  over  time  and  the  initial  opinion 
profile  is  the  one  and  only  fixed  point.  For  a  more  interesting 
scenario,  we  have  the  following  assumption. 

Assumption  8:  The  reward  and  cost  coefficients  satisfy  the 
relation  ^  <  p(0). 

Under  the  strategic  model,  it  follows  from  equation  (12) 
and  Assumption  8  that  the  rate  of  interaction  Pij  changes  over 
time,  depending  on  the  relative  distances  between  agent  i  and 
its  neighbors.  If  the  cost  of  interacting  with  an  agent  exceeds 
the  benefit,  then  the  agent  who  is  to  initiate  an  interaction 
will  impose  a  zero  rate  of  interaction  with  that  agent.  On  the 
contrary,  more  probability  weight  will  be  put  on  the  neighbors 
yielding  higher  (positive)  utilities.  Hence,  the  probability  dis¬ 
tribution  of  pairwise  interactions  Pij  =  {Pij  H-  Pji  )/n  is  also 
dependent  on  the  opinion  distances  between  agents. 

It  follows  from  (9)  and  (12)  that  the  dynamic  of  d  equals 


n 


E 


(15) 


where  r]{d)  =  [ap{d)  —  71]^.  From  (15),  one  can  clearly  see 
that  the  system  stops  evolving  (i.e.,  d  =  0)  if  dij  for  \f{i,j)  G 
E  satisfies  one  of  the  two  conditions:  (i)  dij  =  0;  (ii)  r]{dij)  = 
0.  The  first  condition  is  satisfied  if  the  interacting  agents  are 
in  consensus.  The  second  condition  implies  that  agents  will 
not  interact  if  their  squared  opinion  distance  dij  lies  in  the 
union  Vi  U  [r,  dgup)  where 

=  {de  (0,r)  I  p(d)  <  ^}.  (16) 

Note  that  Vi  is  an  empty  set  when  the  ratio  ji/a  <  p{r~).  In 
this  case,  agents  will  not  update  their  opinions  only  if  dij  >  r 
for  y{i,j)  G  £.  Under  this  situation,  the  network  follows  the 
same  opinion  update  rule  as  in  the  hard-interaction  model  with 
a  threshold  r,  but  it  allows  the  rate  of  interaction  Pij  to  vary 
with  agent  opinion  profiles. 

On  the  other  hand,  as  shown  in  Fig.  3,  when  ^i/a>  p{r~), 
the  set  is  nonempty  and  the  threshold  r  >  inf  (Pi),  where  inf 
denotes  the  infimum.  Since  p{d)  is  a  non-increasing  function 
of  d,  the  range  of  the  set  Pi  goes  from  inf  (Pi)  to  r.  Clearly, 
dij  >  inf  (Pi)  implies  that  the  associated  agents  will  not 
update  their  opinions.  The  opinion  diffusion  in  this  case  will 
not  converge  to  a  consensus  not  only  because  the  agents  may 
not  be  sufficiently  like-minded  (i.e.,  dij  >  r),  but  also  because 
7i  /(T(  is  too  big  to  warrant  sufficient  incentives  between  pairs 
of  like-minded  agents  whose  (squared)  opinion  distances  lie 
in  the  interval  dij  G  [inf(Pi),r). 

2)  Case  Study  (ii):  We  now  examine  the  case  when  the 
energy  function  is  defined  as  ^{d)  =  ^2P^{d)d  for  d  <  r  and 
+(X)  otherwise.  Under  Assumption  2,  it  can  be  shown  that  the 
expression  1/ p{d)  —  l  is  a  non-decreasing  function  of  d,  with 
a  minimum  value  l//i(0)  —  1.  Hence,  if  the  ratio  72/0^  is  less 


Fig.  3.  The  solid  line  represents  the  trust  function  ii{d)  and  the  vertical 
line  maps  the  value  of  fi  such  that  —  p{d)  to  the  point  d  =  inf  (Pi). 

The  “active”  region  means  that  agents  i  and  j  will  interact  if  dij  falls  in  this 
region. 


than  l//i(0)  —  1,  then 


72  ^  I 

ct  p{dij) 


-  1 


for  \/dij.  When  dij  <  r,  replacing  C{dij)  with  ^2p‘^{dij)dij 
in  equation  (12),  observe  that  if  the  preceding  relation  is  true, 
then  Pij  =  0  for  Vj  G  {m  G  A/^  |  dim  <  f}.  This  result, 
together  with  the  fact  that  Pij  =  0  whenever  dij  >  r,  implies 
that  the  initial  opinion  profile  of  the  system  is  the  one  and  only 
fixed  point  because  none  of  the  agents  is  willing  to  interact 
with  its  neighbors.  Hence,  we  make  the  following  assumption. 

Assumption  9:  The  reward  cost  coefficients  satisfy  the  re¬ 
lation  ^  —  1. 

a  —  fi{0) 

We  now  consider  the  opinion  evolution  of  a  system 
satisfying  Assumption  9.  Since  ^{d)  =  j2p‘^{d)d  for 

d  <  r,  the  dynamic  of  the  expected  squared  distance 
d  has  the  same  expression  as  (15)  except  that  p{dij)  = 
\api{dij)  —  (a  +  ^2)p‘^{dij)\  As  mentioned  before,  the  sys¬ 
tem  reaches  a  fixed  point  if  dij  for  V(7,  j)  G  f  satisfies  one 
of  the  two  conditions:  (i)  dij  =  0;  (ii)  p{dij)  =  0.  The  second 
condition  indicates  that  agents  will  not  interact  if  dij  lies  in 
the  union  P2  U  [r,  d^up)  where 

Or  equivalently,  dij  G  (0,  sup(P2)] U [0,  dgup)  since  1  / p{d)  —  1 
is  a  non-decreasing  function  of  d. 

Clearly,  if  the  threshold  r  <  sup(P2)  is  small  relative  to 
the  supremum  of  the  set  V2,  then  agents  in  the  network  will 
not  update  their  opinions  because  either  they  are  too  closed- 
minded  to  the  opinions  of  the  others  or  they  do  not  have 
sufficient  incentives  to  interact.  In  contrast,  as  depicted  in 
Fig.  4,  if  the  threshold  r  sup(P2),  agents  will  interact 
if  dij  lies  in  the  open  interval  (sup(P2),  f).  Hence,  it  can  be 
deduced  that  the  system  will  form  one  or  multiple  opinion 
clusters.  Within  each  cluster,  the  squared  opinion  distances 
are  upper  bounded  by  sup(P2)-  Between  the  clusters,  the 
squared  opinion  distances  are  lower  bounded  by  the  threshold 
r.  Comparing  Figures  3  and  4,  we  observe  that  the  two  cost 
functions  lead  to  distinct  active  regions  (under  appropriate 
conditions  on  the  cost-reward  coefficients):  one  case  can  lead 
to  consensus;  the  other  case  leads  to  clustering. 
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Fig.  4.  The  solid  line  represents  the  trust  function  and  the  vertical  line 
maps  the  value  of  ^  such  that  72/0;  =  1/ —  l  to  the  point  d  =  sup(D2)- 

V.  Numerical  Results 

In  this  section,  we  first  numerically  validate  the  analytical 
results  of  both  the  hard-interaction  model  (Sec.  V-A)  and 
the  strategic  interaction  model  (Sec.  V-B)  using  synthetic 
data.  Then  we  offer  evidence  that  the  proposed  model  for 
opinion  evolution  could  portray  real  phenomena  using  the 
Social  Evolution  dataset  [23]  (see  Sec.  V-C). 

A.  Hard-Interaction  Model 

Recall  that  the  choice  of  the  underlying  communication 
network  Q  is  arbitrary  and  the  analytical  results  hold  for  any 
connected  network.  We  start  by  generating  Q  using  a  random 
geometric  graph  (RGG),  i.e.,  Q  =  Q{n,r),  which  consists 
of  n  randomly  distributed  social  agents  over  an  unit  disk 
with  a  radius  of  communication  r.  While  r  =  1  implies 
a  fully  connected  network,  r  =  0  represents  a  completely 
disconnected  network.  Each  initial  opinion  profile  x^O]  for 
Vi  G  V  is  uniformly  distributed  in  the  opinion  space  iT  with 
q  =  3  possible  decision  states.  Without  loss  of  generality, 
the  ^2 -norm  is  used  to  measure  the  opinion  distance  between 
agents,  i.e.,  s(x^,Xj)  =  ||xi  —  Xj||^  with  A  =  Iq. 

Eor  the  hard-interaction  model,  we  set  the  trust  function  to 
be  /i(d)  =  0.4  for  d  <  r  and  0  otherwise.  Three  connected 
RGG  graphs  0(n,O.8)  are  generated  with  n  =  50,100,200 
number  of  agents  and  are  kept  fixed  for  all  trials.  We  define 
a  subgraph  Gcfflk]  =  (V,feff[^])  of  G  at  each  time  instant, 
where  feff[^]  =  {{hj)  ^  ^  I  dij[k]  <  r}.  A  normalized 
algebraic  connectivity  of  the  subgraph  Gcffik]  is  defined  as 
the  algebraic  connectivity  of  Gcffik]  divided  by  the  algebraic 
connectivity  of  G^  Eor  k  sufficiently  large,  the  normalized 
algebraic  connectivity  approaches  either  0  or  1.  A  value  of  0 
indicates  that  ^eff  is  disconnected;  A  value  of  1  indicates  that 
^eff[^]  =  ^  and  moreover,  dij[k]  ^  0  for  V(7,  j)  G  f ,  implying 
that  consensus  is  reached.  The  top  panel  in  Eig.  5  shows  the 
(normalized)  algebraic  connectivity  of  the  subgraph  Geff[k], 
where  k  is  sufficiently  large,  averaged  over  300  realizations  for 
various  values  of  r.  Each  realization  starts  with  an  uniformly 
distributed  initial  opinion  profile.  Uniform  communication  rate 
(i.e..  Pi  =  1/n  and  Pij  =  l/{n  —  1))  is  used.  Observe  that 
as  r  approaches  0.64,  society  reaches  a  consensus  almost 
surely  and  thus  a  phase  transition  occurs  approximately  at 
r  =  0.64.  Moreover,  as  n  increases,  the  phase  transition 
boundary  becomes  sharper;  the  normalized  algebraic  connec¬ 
tivity  decreases  at  r  =  d[0].  The  bottom  panel  in  Eig.  5 
shows  the  corresponding  histograms  of  the  300  asymptotic 


>. 
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Fig.  5.  Phase  transition  (top)  of  a  society  from  polarized  beliefs  to  a 
consistent  belief.  Histogram  (bottom)  of  300  asymptotic  opinion  profiles  with 
n  =  100  at  r  =  0.16,  0.36,  0.64. 

opinion  profiles  at  r  =  0.16,0.36,0.64,  respectively.  Notice 
the  evidence  of  fragmentation  or  clustering  at  smaller  values 
of  r,  and  consensus  at  larger  values. 

Eig.  6  compares  the  phase  transitions  for  three  different 
network  configurations:  RGG,  Erdos-Renyi  graph  and  small- 
world  graph.  Specifically,  n  =  50  agents  are  generated  for 
each  of  the  graphs.  The  radius  of  communication  r  in  the 
RGG  is  set  to  be  0.8.  The  Erdos-Renyi  graph  [25]  is  generated 
uniformly  at  random  from  the  collection  of  all  graphs  which 
have  n  =  50  nodes  and  M  =  120  edges.  The  small- 
world  graph  [26]  is  generated  by  setting  psw  =  1,  Qsw  =  1 
and  Gsw  =  2,  where  Qsw  denotes  the  maximum  distance 
of  short-range  connections,  Qs^  represents  the  number  of 
random  connections  to  add  per  node  and  Gqw  is  the  clustering 
exponent.  All  of  the  three  graphs  used  in  the  simulation  are 
connected.  Uniform  rate  of  communication  Pij  is  adopted. 
The  (normalized)  algebraic  connectivity  of  Gcffik]  for  each 
graph  is  averaged  over  400  realizations,  each  of  which  starts 
with  an  uniformly  distributed  initial  opinion  profile.  Observe 
from  Eig.  6  that  the  phase  transitions  almost  overlap  with  one 
anther.  As  discussed  towards  the  end  of  Section  III,  specific 
choices  of  the  fixed  graph  (P)  do  not  impact  the  dynamics 

B.  Strategic  Interaction  Model 

Eor  the  strategic  interaction  model,  we  set  /^  =  2  and 
define  the  trust  function  to  be  11(d)  =  0.5  —  ^.Ad?  for 
Vd  <  r  as  shown  in  Eig.  7  for  different  values  of  r,  i.e., 
r  =  0.25,  0.64,  0.81.  The  underlying  communication  graph  is 
a  RGG  graph  with  n  =  50  and  r  =  0.8.  The  initial  opinion 
profiles  for  the  individual  agents  are  generated  uniformly  over 
the  opinion  space  A’  with  q  =  3.  Opinion  distances  between 
agents  are  measured  by  the  ^2 -norm. 

I)  Case  Study  (i):  Consider  the  subgraph  GeE[k]  = 
(V,feff[^])  of  G  at  time  k,  where  contains  all  the 

edges  (i,j)  G  E  whose  corresponding  distances  dij[k]  are 
such  that  dij[k]  <  r  if  Pi  =  0  and  dij  <  inf  (Pi)  otherwise. 
Eig.  8  shows  the  (normalized)  algebraic  connectivity  of  the 
graph  GeE  [k]  for  k  sufficiently  large.  The  plot  is  averaged  over 
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Fig.  6.  Phase  transitions  using  three  different  random  graphs,  n  =  50. 


Fig.  7.  Trust  Function  with  r  =  0.25,  0.64,  0.81 


400  realizations  for  different  values  of  r  and  inf  (Pi).  Each 
realization  starts  with  an  uniformly  distributed  initial  opinion 
profile.  Observe  that  when  inf  (Pi)  is  small,  i.e.,  71/ct  is  large, 
the  agents  are  less  likely  to  reach  a  consensus  for  any  value 
of  r.  In  contrast,  when  inf  (Pi)  is  large,  i.e.,  inf  (Pi)  >  0.64 
approximately,  the  society  tends  to  form  a  convergent  opinion 
almost  surely  for  large  values  of  r  (approximately  above  0.64). 

2)  Case  Study  (ii):  Fig.  9  shows  the  final  outcome  of  the 
interactions  (top  panel)  and  the  squared  distance  distribution 
(bottom  panel)  with  r  =  0.09  and  sup(P2)  =  0.0158  (i.e., 
72/0^  =  4.0016).  Observe  from  the  top  panel  that  three 
opinion  clusters  are  formed.  Within  each  cluster,  the  (squared) 
opinion  distances  are  upper  bounded  by  sup(P2)  =  0.0158, 
as  shown  in  the  bottom  right  panel  of  Fig.  9.  Also,  the 
(squared)  distances  between  clusters  are  at  least  0.18,  which 
is  much  larger  than  the  threshold  r  =  0.09.  Fig.  10  shows  the 
final  outcome  of  the  interactions  (top  panel)  and  the  squared 
distance  distribution  (bottom  panel)  when  r  =  0.64  and 
sup(P2)  =  0.0158.  In  this  case,  agents  form  a  single  opinion 
cluster,  as  shown  in  the  top  panel;  the  squared  distances  within 
this  cluster  are  computed  to  be  upper  bounded  by  0.0156, 
which  is  less  than  sup(P2)-  Fig-  H  shows  the  time  evolution 
of  the  expected  squared  opinion  distance  d[k]  as  k  increases 
for  the  two  scenarios.  In  both  cases,  d  decreases  until  it  reaches 
an  equilibrium  point. 

C.  Social  Evolution  Dataset 

To  demonstrate  the  validity  of  the  proposed  interaction 
models,  we  use  data  collected  in  the  Social  Evolution  ex¬ 
periment  [23].  Specifically,  the  experiment  monitored  more 
than  80%  of  the  residents  in  a  dormitory,  with  a  population 
of  approximately  30  freshmen,  20  sophomores,  10  juniors,  10 
seniors  and  10  graduate  student  tutors.  Interactions  between 
individuals  were  tracked  by  their  proximity,  location,  SMS 
and  call  records.  Surveys  were  conducted  monthly  on  so¬ 
cial  relationships,  health-related  habits,  on-campus  activities, 
political  views  and  common  cold  symptoms.  The  friendship 
or  communication  network  is  connected;  the  exact  network 
graph,  although  not  significant  to  this  research,  is  provided  in 


Fig.  8.  Phase  Transitions  for  different  values  of  inf(Di) 
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Fig.  9.  Final  opinion  landscape  (top),  initial  (bottom  left)  and  final  (bottom 
right)  squared  distance  distributions,  with  r  =  0.09  and  sup(I72)  =  0.0158. 


[27]  and  we  assume  it  stays  constant.  While  many  insights  can 
be  gleaned  from  an  analysis  of  the  full  data,  for  this  research, 
the  data  that  we  choose  to  exploit  are  the  location  data  and 
the  survey  data  on  health-related  habits. 

The  location  data  are  collected  through  a  mobile  phone 
application  that  scans  the  wireless  local  area  network  (WLAN) 
access  points  within  a  certain  range.  We  use  the  location  data 
to  study  how  social  interactions  shape  individuals  preferences 
in  terms  of  how  frequently  they  visit  a  place.  There  are  32724 
identified  WLAN  access  points;  we  choose  to  use  the  top 
100  most  frequently  visited  access  points  (approximately  10 
locations)  to  extract  individual  opinions.  The  opinion  space  X 
is  of  dimension  =  101,  where  the  first  100  states  correspond 
to  the  top  100  frequently  visited  WLAN  access  points  and  the 
last  state  represents  all  other  locations.  Opinions  are  estimated 
once  every  four  weeks.  Specifically,  we  assume  a  Binomial 
distribution  for  visiting  an  access  point  with  being  the 
total  number  of  places  a  student  could  visit  in  four  weeks. 
We  count  the  total  number  of  times  a  student  visited  one 
of  the  100  access  points.  Then  the  probability  of  visiting 
each  place  is  the  maximum  likelihood  estimator  (MSE)  of 
the  Binomial  distribution,  given  the  actual  count.  The  relative 
number  of  visits  to  different  location  is  an  indication  of  the 
importance  of  the  locations.  A  question  of  interest  is  how 
this  relative  importance  evolves  over  time,  and  whether  this 
stabilizes.  Opinion  distances  are  measured  by  the  ^2 -norm 
with  A  =  Iq.  Although  the  distance  calculations  are  carried 
out  with  q  =  101,  it  is  important  to  note  that  the  intrinsic 
dimensionality  of  the  opinion  space  A'  can  be  estimated  by 
nioc  H-  1  where  nioc  is  the  number  of  locations  with  which 
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Fig.  10.  Final  opinion  landscape  (top),  initial  (bottom  left)  and  final  (bottom 
right)  squared  distance  distributions,  with  r  =  0.64  and  sup(P2)  =  0.0158. 
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Fig.  11.  Dynamics  of  expected  (squared)  distance  d  for  r  =  0.09  (left)  and 
r  =  0.64  (right). 

the  WLAN  access  points  are  associated;  according  to  [27], 
nioc  ~  10  for  these  top  100  WLAN  access  points. 

Fig.  12  shows  the  histograms  of  the  squared  opinion  dis¬ 
tances  d  in  the  following  six  time  periods:  (PI)  10/1/08  - 
10/25/08;  (P2)  10/26/08  -  11/22/08;  (P3)  11/23/08-12/20/08; 
(P4)  12/21/08-1/17/09;  (P5)  1/18/09-2/14/09;  (P6)  2/15/09- 
2/28/09.  We  observe  that  in  PI,  d  can  take  almost  all  the  values 
ranging  from  0  to  about  0.7  and  a  few  values  around  1,  while 
by  P6,  the  values  of  d  are  concentrated  around  0  and  1.  We 
compute  a  quantity  |cc|o,  the  number  of  “quantized”  distances 
d  that  have  nonzero  probability:  the  value  Iccjo  drops  from 
59  in  PI  to  17  in  P6,  indicating  that  opinions  are  clustered. 
Fig.  13  shows  the  time  evolution  of  the  expected  (squared) 
opinion  distance  d.  We  observe  a  similar  downward  trend  in 
d  (see  Fig.  11)  except  in  the  weeks  around  Christmas,  i.e.,  P4. 
The  dotted  line  in  Fig.  13  shows  a  predicted  curve  for  d  once 
the  data  in  P4  were  discarded. 

Consider  now  the  monthly  survey  data  on  health-related 
habits,  in  particular,  the  eating  habits  and  the  exercising  habits. 
Eating  habits  are  tracked  by  the  number  of  healthy  items  a 
person  eats  (i.e.,  number  of  salads  per  week  and  number  of 
fruits  and  vegetables  per  day).  Exercising  habits  are  recorded 
by  the  number  of  aerobic  exercises  (each  lasting  at  least  20 
minutes)  per  week  and  the  number  of  times  one  participates 
in  team  sports  per  week.  Thus  individual  opinions  can  be 
constructed  by  computing  the  probabilities  of  the  following 
four  states:  (1)  eating  healthy  but  not  exercising;  (2)  exercising 
but  not  eating  healthy;  (3)  eating  healthy  and  exercising;  (4) 
none  of  the  above.  The  opinions  are  again  computed  on 
a  monthly  basis.  Assume  that  both  events  (eating  healthy 
and  exercising)  are  independent  and  each  follows  a  Binomial 
distribution  with  ni  and  n2  being  the  total  number  of  items  a 
person  could  eat  per  month  and  the  total  number  of  exercise 


Squared  Opinion  Distance  d  Squared  Opinion  Distance  d 


Fig.  12.  Histograms  of  squared  opinion  distance  d  in  PI  and  P6.  The  Iq 
norm  |a;|o  denotes  the  number  of  “quantized”  opinion  distances  d  that  have 
nonzero  probability. 


Fig.  13.  Evolutions  of  the  expected  d:  measured  v.s.  predicted  dynamics 
from  PI  to  P6. 

activities  one  could  do  per  month.  Then  the  probability  of 
each  event  is  the  maximum  likelihood  estimator  (MSB)  of 
the  Binomial  distribution  given  the  actual  number  of  healthy 
items  a  person  eaten  in  a  month  or  the  number  of  exercises  a 
person  done  in  a  month.  Using  the  independence  assumption, 
the  probabilities  of  the  four  states  can  be  computed. 

Fig.  14  shows  the  histograms  of  the  squared  opinion  dis¬ 
tances  in  selected  months'^.  Clearly,  they  show  no  significant 
change  over  time.  Unlike  the  previous  experiment  (on  the 
frequently  visited  places)  which  exhibits  clustering  behavior 
as  time  evolves,  social  interactions  between  students  do  not 
seem  to  affect  individual  health-related  habits.  This  suggests 
that  social  interactions  do  not  influence  a  person’s  opinions 
toward  different  issues  at  the  same  level,  in  terms  of  the  rate 
of  behavioral  changes.  There  are  two  possible  explanations. 
One  is  that,  generally  speaking,  health-related  habits  usually 
change  over  a  long  time  scale  and  this  might  be  caused 
by  a  small  threshold  r  in  the  trust  function;  people  are 
close-minded  toward  changing  habits.  Hence,  it  is  difficult 
for  the  survey  data  to  capture  such  minuscule  changes  in 
opinion.  Another  reason  is  that  students  are  often  not  actively 
involved  in  exchanging  relevant  health  advices  during  social 
interactions.  If  this  is  the  case,  then  Pij  is  in  fact,  equal  to 
zero  for  i  ^  j  and  thus  opinions  stay  constant. 

VI.  Conclusions 

We  proposed  a  generalization  of  the  Deffuant-Weisbuch 
model  and  studied  opinion  dynamics  in  a  connected  network 
under  two  related  interaction  models,  i.e.,  the  hard-interaction 
model  and  the  strategic  interaction  model.  Under  the  hard- 
interaction  model,  we  provided  a  necessary  condition  that 
guarantees  opinion  convergence.  We  showed  that  the  necessary 
condition  for  convergence  does  not  change  as  long  as  the 

^Survey  results  are  not  available  for  Nov.  08,  Jan.  09,  Feb.  09  and  May  09. 
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Fig.  14.  Histogram  of  squared  opinion  distance  d  for  selected  months. 


communication  rates  are  time-invariant.  Under  the  strategic  in¬ 
teraction  model,  we  showed  how  opinion  formation  processes 
are  affected  by  the  individual  incentives  behind  interactions. 
In  particular,  we  explored  two  specific  utility  functions,  which 
lead  to  two  different  asymptotic  opinion  patterns.  Our  analysis 
of  a  rich  social  data  indicate  evolution  towards  clustering 
behaviors,  measured  by  the  frequency  of  visiting  a  location; 
however,  the  data  do  not  indicate  clustering  with  respect  to 
other  behavioral  traits  which  may  evolve  at  a  longer  time  scale. 

VII.  Appendix 

Derivation  of  Eqn.  (8):  Given  that  agents  i  and  j  interacted 
at  /c  -t-  1,  the  change  in  D  after  the  interaction  is  only  affected 
by  the  opinion  distances  between  the  pairs  (i,  j),  {i,£)  and 
(j,  i)  for  V^.  It  then  follows  from  (5),  (6)  and  (7)  that  [k  -|- 
1]  —  D[k]  equals 

dij  “h  1]  dij 

+  diilk  +  1]  —  dii[k]  -f  dji[k  +  1]  —  dji[k] 

=  -  4fi{dij[k]){l  -  iJ.{dij[k]))dij[k] 

+  [2ij,‘^{dij[k])dij[k] -2ij,{dij[k])dij[k]] 

=  -  4:ij,{dij[k]){l  -  id{dij[k]))dij[k] 

+  (n  -  2)  [2 11^ {dij  [k]) dij  [k]  -  2ij,{dij[k])dij[k]] 

=  -  2ni2{dij[k])  [1  -  fi{dij[k])]  dij[k]. 

Lemma  5:  Let  h{d)  =  p{d)d  =  /i(d)(l  —  p{d))d.  Then 
under  Assumptions  2  and  4,  its  second  derivative  equals  h  = 
— 2(/i)^(i+(l  — — 2/i)/i  <  0  for  d  G  (0,r)  and  thus 
h{d)  is  concave  for  d  G  (0,r).  Let  Peff  =  ft {u)du  <  1 
and  deff  =  uft{u)du.  Then 

E{h{d)}  =  [  h{u)ft{u)du  =  Pefi  [  h{u)^^^du. 

Jo  Jo  ^eff 

Using  Jensen’s  inequality  and  the  relation  deir  <  d  yields 
E{/i(<i)}  <  Pesh  <  P  ^  -  pi^es)d. 

The  last  inequality  holds  because  p{d)  is  a  non-increasing 
function  of  d  (i.e.,  its  derivative  p  =  {1  —  2p)fi  <  0).  We 


Fig.  15.  Vector  Field;  the  point  ^(s;  t)  moves  toward  the  stable  fixed  point 
G  =  0  if  ^(0;  t)  <  K  and  moves  away,  otherwise. 


have  p{deg)/p{d)  <  p{0)/p{T  )  <  /3  and  subsequently, 
E{/i(d)}  <  I3h{d). 

Lemma  6:  For  s  sufficiently  small,  Taylor’s  expansion 
gives  p  {b{t  +  s))  K  p  {b{t))  +  {b{t  +  s)  -  b{t))  p  {b{t)).  Then 
system  b{t  +  s)  =  —I3p{b{t  +  s))b{t  +  s)  becomes 


b{t  +  s) 


-(3p{b{t))b{t  +  s)  if  p  (6(f))  =  0; 
-f{s-,t)  if  p  (6(f))  <0. 


where  /(s;f)  =  /3  p  {b{t))  -  b{t) p  {h{t))  (l  -  6(f  + 

s)  and  p{h)  :=  dp/dh.  In  the  first  case,  p{b{t))  is  locally 
constant  and  hence,  the  local  rate  of  convergence  around  b{t) 
is  exponential  and  is  equal  to  p{b{t))  when  p{b{t))  >  0, 
i.e.,  b{t)  <  r.  For  the  second  case  when  p{b{t))  <  0,  define 
^(5;  t)  :=  b{t  +  s)/b{t).  With  respect  to  s,  the  dynamics  of  ^ 
become 


C(s;  t)  =  +  s)  =  -l3R{t)^{s;  t) 


where  R{t)  =  p{b{t))  —  b{t)p{b{t)) 

Pirn) 


m) 


K{t)  =  1  + 


-b{t)p{b{t)) 


>  1 


(18) 

(19) 

(20) 


Note  that  the  dynamics  of  resemble  the  logistic  equa¬ 

tion.  There  are  two  equilibria  at  =0  (stable)  and  ^2  =  K {t) 
(unstable).  When  R{t)  >  0,  as  shown  in  Fig.  15,  the  system 
will  converge  if  K{t)  >  <^(0;t)  =  1.  Indeed,  if  b{t)  <  r, 
then  K{t)  =  1  —  ^  converges.  The 

exponential  rate  of  convergence  equals  R{t).  Moreover,  since 
6(t)  is  a  monotonically  decreasing  function  of  t,  an  equivalent 
condition  for  convergence  is  6(0)  <  r.  On  the  contrary,  when 
b{t)  =  r,  we  get  K{t)  =  1  because  p(r)  =  0.  For  b{t)  >  r, 
both  p  {b{t))  and  p  {b{t))  are  zero,  which  implies  R{t)  =  0  in 
(18)  and  bifurcation  occurs.  Thus  the  system  may  not  converge 
when  b{t)  >  r  for  Vt. 
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